Random Fields and the Partially-Paramagnetic State of CsCoo.83Mgo.i7Br 3 : 

A Critical Scattering Study 
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Critical neutron scattering measurements have been performed on CsCoo.83Mgo.i7Br3, a dilute 
stacked triangular lattice (STL) Ising antiferromagnet (AF). A two component lineshape associ- 
ated with the critical fluctuations appears at a temperature coincident with Tati observed in pure 
CsCoBr3. Such scattering is indicative of fluctuations in prototypical random field Ising model 
(RFIM) systems. The random field domain state arises in this case due to geometrical frustration 
within the STL Ising AF, which gives rise to a 3 sublattice Neel state, in which one sublattice is 
disordered. Magnetic vacancies nucleate AF domains in which the vacancies reside on the disordered 
sublattice thereby generating a RFIM state in the absence of an applied magnetic field. 
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Geometrically frustrated magnetic materials have been 
of intense recent interest, in part due to the diversity in 
the exotic ground states that they display Q. These 
span the range from complex non-collinear Neel states, 
to disordered co-operative paramagnetic and spin liq- 
uid states, to spin glass-like states in the presence of 
little or no chemical disorder. Among the exotic mag- 
netic states known to exist are partially-paramagnetic 
Neel states, characterized as multi-sublattice antiferro- 
magnets, in which one of the sublattices in this long range 
ordered structure remains paramagnetic. This unusual 
structure occurs in several magnetic materials made up 
of triangular planes, including the stacked triangular lat- 
tice Ising antiferromagnetic insulators CsCoBr3 || and 
CSC0CI3 H , as well as the hexagonal intermetallic com- 
pound UNi 4 B {§. 

A striking feature of the phase behavior of CsCoBi'3 
and CSC0CI3 is their extreme sensitivity to the presence 
of non-magnetic impurities on the magnetic Co site. Neu- 
tron scattering work on CsCoi-^Mg^Cls by Nagler et 
al. H and on CsCoi-^Mg^Brs by Rogge et al. || show 
that 17% dilution of the cobalt sites with non-magnetic 
magnesium appears to wipe out at least two (and possi- 
bily three) long range ordered states which pure CSC0CI3 
displays below T^i^ 21 K ||, and which pure CsCoBr 3 
displays below ~ 28 K §. This is surprising as 3 dimen- 
sional lattices support percolating long range order to ~ 
75% dilution. §. 

In this letter, we report on a critical neutron scatter- 
ing study of single crystal CsCoo.83Mg .i7Br 3 . We show 
that the critical fluctuations contain two components, the 
narrower of which (in Q-space) displays an onset at a 
temperature co-incident with Tjvi of pure CsCoBr3, ~ 
28 K. Two distinct components to the fluctuations are 
easily observed over the approximate temperature inter- 
val for which pure CsCoBr3 exhibits the exotic partially- 
paramagnetic, 3-sublattice Neel state. We argue that this 
two component lineshape is characteristic of a domain 
state typical of the Random Field Ising Model (RFIM), 



and which is generated by the presence of quenched mag- 
netic vacancies and the paramagnetic sublattice. 

RFIMs occupy interesting and important territory 
within the physics of disordered media ||. Within this 
model the disorder arises due to random local symmetry 
breaking, as opposed to random interactions and frus- 
tration, which drives the physics of conventional spin 
glasses M. The RFIM is known to describe well the 
properties of magnetically-dilute Ising antiferromagnets 
in the presence of magnetic fields jl(| O, |l2| [l3] , as well 
as a wide class of other disordered materials, such as, 
for example, phase separating fluids in the presence of 
gels @. 

Imry and Ma first considered the effects of a ferro- 
magnetic Ising system in the presence of a site-random 
magnetic field, and argued that the random fields would 
destroy the long range order of such a classical system for 
dimensions less than 4 ]l5| . Fishman and Aharony later 
argued that magnetically dilute Ising antiferromagnets 
in the presence of a uniform magnetic field were isomor- 
phous to the Imry-Ma scenario |l6| . This was important 
as not only do site-random antiferromagnets occur in na- 
ture, but it is easy to apply an external magnetic field 
and thereby tune the amplitude of the random field. 

Systematic studies of such site-random Ising antiferro- 
magnets have been carried out for both three dimensional 
materials (such as CoZni_ x F 2 friof , MnZni_ x F2 juj, and 
FeZni_ ;E F2 1 12 ) as well as quasi-two dimensional sys- 
tems (such as Rb2Coo.7Mgo.3F4 1 13 ). While some open 
questions remain, it is clear that the field-cooled state 
at low temperatures corresponds to a system broken up 
into many relatively small domains. A characteristic 
signature of such a state in neutron scattering experi- 
ments is a two-component lineshape in wave vector, with 
a sharp and a broad component, as opposed to the single- 
component lineshape corresponding to, say, Ornstein- 
Zernike fluctuations [l7| . Such a two-component line- 
shape, however, was inferred from the results of fitting 
the neutron lineshapes to an appropriate addition of two 
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FIG. 1: Role of non-magnetic atoms in generating a random 
field domain state in CsCoo.83Mgo.i7Br3 



terms. 

CsCoo.83Mgo.i7Br3 belongs to the ABX3 family of anti- 
ferromagnetic (AF) insulators |Q . The interest in these 
materials stems from the quasi-one-dimensional nature 
of their magnetic interactions, as well as their geometri- 
cally frustrated cooperative behaviour. Pure CsCoBr 3 is 
an example of an Ising-like antiferromagnet on a stacked 
triangular lattice gj. It crystallizes in the hexagonal 
space group PG^/mmc, with low temperature lattice pa- 
rameters a= 7.45 and c= 6.30 A. The Co 2+ ions possess 
spin-i magnetic moments, which are aligned along the 
hexagonal c direction, and lie on a simple hexagonal lat- 
tice. Inelastic neutron scattering studies have shown that 
this material is described by the following spin Hamilto- 
nian @: 

H = ^(2J{5f^ 1 + e(5f5f +1 + ^ +1 )} 



+h s!{~ir + 2j'Y J s z l s z J ) 



(1) 



where J= 1.62 THz, J'= 0.0096 THz and e= 0.137. 
The dominant AF interactions along the c-axis result in 
quasi-one dimensional behaviour over an extended tem- 
perature regime. At sufficiently low temperatures, how- 
ever, the weaker ab-plane interactions cause two or pos- 
sibly three, magnetic phase transitions to three dimen- 
sional long range ordered Neel states ^ pfj| . The first 
phase transition, from the paramagnetic state to a long 
range ordered, but partially-paramagnetic state, occurs 
at Tati=28.3 K. This state is characterized by three sub- 
lattices, in which two out of every three spins on trian- 
gular plaquettes are ordered antiferromagnetically, while 
the third remains disordered. This is shown as +, - and 
P in Fig |} The ordering triples the periodicity within 
the ab-plane, giving rise to magnetic Bragg reflections of 
the form (-|, -|, I) with 1 odd. Below ~ 13 K the param- 
agnetic sublattice orders fully, either (+) up or down (-), 
to form a ferrimagnetic sheet within the ab-plane. The 
net moment within this ferrimagnetic sheet orders alter- 
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FIG. 2: Peak intensity of the magnetic (§, §, 1) and (|, |, 1) 
Bragg reflection of CsCoBr3 and CsCoo.83Mgo.i7Br3 respec- 
tively as a function of temperature. 



nately from ab-plane to ab-plane, so as to produce an 
antifcrromagnetic structure, with no net magnetic mo- 
ment. An intermediate magnetically-ordered structure 
may exist between ~ 18 K and 13 K, in which the para- 
magnetic sublattice gradually orders. 

Classical Ising antiferomagnets on stacked triangu- 
lar lattices have been studied extensively by Landau- 
Ginsburg- Wilson 21 and Monte Carlo techniques [p2| . 
The high and low temperature long range ordered states 
are accounted for by all theoretical models. 

The argument by which the combination of a partially- 
paramagnetic Neel state and magnetic vacancies gener- 
ate a random field domain state is illustrated in Fig. 
1. These quenched vacancies will preferentially nucle- 
ate domains in which the vacancies are coincident with 
the paramagnetic sublattice so as to minimize the loss 
of exchange energy. Two such domains are depicted in 
Fig. 1, with a domain wall running between them. In 
the absence of the paramagnetic sublattice, the vacan- 
cies are co-incident with either up or down sublattices, 
and only break symmetry between the two in the pres- 
ence of an externally applied longitudinal magnetic field; 
the Fishman-Aharony scenario ||16fl . 

The critical neutron scattering study was carried out 
on the N5 triple axis spectrometer, operated in two axis 
mode, at the NRU reactor of Chalk River Laboratories. 
The high quality single crystal was the one studied by 
Rogge et al. Q. Neutrons of wavelength 2.37 A were 
reflected from the (0, 0, 2) planes of a pyrolitic graphite 
(PC) monochromator. A PC filter was placed between 
the monochromator and sample to remove higher order 



contamination. The sample, mounted in a He cryostat, 
was oriented with its (h, h, I) plane in the horizontal scat- 
tering plane. Collimation of 0.19° in the incident beam 
and 0.25° in the diffracted beam was used to achieve rel- 
atively high Q resolution. 

The temperature dependence of the Bragg 
peak intensities in both pure CsCoBr3 and dilute 
CsCoo.83Mg .i7Br 3 is shown in Fig. |[ The square of the 



order parameters were taken to be proportional to the 
peak intensities at the magnetic Bragg reflection, (|, |, 1) 
for CsCoBr 3 and (|, |, 1) for CsCoo.83Mgo.i7Br 3 . 

While more complete than that reported by Rogge et 
al. [|], the CsCoo.83Mgo.i7Br3 order parameter squared 
presented in Fig. 2 reproduces the key characteristics of 
the earlier measurements. The Bragg intensity displays 
upwards curvature down to at least 5 K, and the clear 
phase transitions in pure CsCoBr3 at ~ 28.3 K and ~ 13 
K are not evident. Mean field behavior would predict this 
intensity to grow linearly with decreasing temperature as 
~ (Tp-T) 2 ^ with /3=0.5. Typical three dimensional crit- 
ical behavior would give /3~ |. Pronounced downwards 
curvature in the Bragg intensity as a function of temper- 
ature would be expected, and indeed is observed in pure 
CsCoBr3 as shown in Fig. 2. 

Similar order parameter measurements on 
CsCoo.83Mgo.17d3 by Nagler et al. || also show 
upwards curvature to the Bragg intensity as a function 
of temperature to temperatures below 5 K, and again no 
evidence of the Tjvi^ 21 K observed in pure CSC0CI3. 
While these earlier studies concluded phase transitions 
occur near 9 K and 6 K, for CsCoo.83Mgo.i7Br3 and 
CsCoo.83Mgo.17d3 respectively, neither study showed 
compelling evidence for any phase transition above ~ 5 
K. 

As a function of temperature, scans in h and I were 
performed through three ordering wavevectors: Q rd= 
(|, |, 1),(§, |, 1) and (|, |, 1). As shown by the h-scans 
in Fig. [| a sharp, albeit relatively weak component, ad- 
ditional to the broad critical fluctuations, appears below 
~ 28 K. This indicates that a broken symmetry has oc- 
curred very close to pure CsCoBr3's Tjvi- Over the tem- 
perature range ~ 28 K to ~ 20 K, the scattering clearly 
displays a two component line shape. 

The data were analysed by convoluting the scattering 
cross section with the measured 3-dimensional resolution 
function. The cross section, centered at each ordering 
wavevector, is described by the sum of a broad and a nar- 
row Lorentzian, both having anisotropic widths. The cal- 
culated intensities included both the Co 2+ magnetic form 
factor, and the neutron polarization factor, by which only 
components of moment perpendicular to Q contribute to 
the scattering. 

The inverse correlation lengths within the ab-plane 
and along the c-axis were allowed to differ, acounting 
for the strong anisotropic nature of the AF interactions 
within this material: hence: 
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where < n' ab , k c < k' c , and q = Q — Q or d- The inverse 
correlation lengths of the broad Lorentzian were allowed 
to refine, while those of the narrow one were held fixed (at 
K a b=0- 00546 and k c =0. 00624 A -1 , approximating their 
resolution limit). A complete fit to the data required 
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FIG. 3: (top) Scattering data along the (h, h, 1) direction at 
the 13 measured temperatures. The white arrows indicate the 
temperatures at which clear phase transitions are observed in 
pure CsCoBr3 (see Fig. 2) . (bottom) Representative (h, h, 1) 
scans clearly showing the onset of a two-component line shape 
below Tjvi. The solid lines are fits of the data to Eq. 2, 
convoluted with the experimental resolution. The bottom 
data set, scaled in intensity by a factor of 10 3 , shows that 
only the sharp component remains at 5 K. 



adding two broad and relatively weak Lorentzians cen- 
tered at wavevectors (|,|,1) and (|,|,1), in addition 
to a constant background. This extra scattering origi- 
nates from the domain walls within the ab-plane. A sta- 
ble paramagnetic sublattice requires a net zero exchange 
field from its near neighbours (3 up (+) and 3 down (-) 
near neighbours). The domain wall regions therefore ex- 
clude the paramagnetic sublattice and reflect a local two 
sublattice structure (Fig. [l]), giving rise to scattering at 
(§,§,/) with I odd. 

Figure || shows the temperature dependence of the re- 
fined parameters. The narrow component to the line- 
shape turns on below 28 K, while the broad component 
extents above Tjvi- Below 20 K the intensity of the nar- 
row component rapidly increases, and the data can be 
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FIG. 4: (top) Temperature dependence of the amplitude of 
both the narrow and broad Lorentzians, resulting from fits 
shown in the bottom panel of Fig. 3. (Bottom), The tem- 
perature dependence of the inverse correlation lengths of the 
broad Lorentzian are shown along with fits to divergences in 
the correlation lengths, appropriate to both mean field theory 
(^=i) and the 3D Heisenberg model (i/=0.707), as described 
by Eq. 3 



equally well described by a single component Lorentzian 
with anisotropic line widths. 

To seek out a second, lower phase transition, we have 
fit the temperature dependence of the inverse correlation 
lengths of the broad component between 20 K and 28 
K. This was fit to a standard critical divergence of the 
correlation lengths as: 



(t - T N2 y 



(3) 



We fit with critical exponents v ranging from a low of A , 
appropriate to mean field theory, to a high of 0.707 ap- 
propriate to the 3D Heisenberg model. These fits, shown 
in the bottom panel of Fig. 4, are consistent with a sec- 
ond phase transition occuring between ~ 16 K and 18 
K. 

While a unique, quantitative analysis of the scatter- 
ing below 20 K is difficult, examination of the qualitative 
features of the scattering in Fig. 3 is useful. The order 
parameter for CsCo .83Mgo.i7Br 3 in Fig. 2 shows little or 
no indication of the ~ 13 K transition in CsCoBr'3. How- 
ever the (h,h,l) scans in the top panel of Fig. 3 clearly 
show a narrowing in the diffuse scattering below 13 K, 



indicating the onset of a true long range ordered state. 
We interpret this as a tendency towards full ordering of 
the paramagnetic sublattice below T7V3 ~ 13 K, thereby 
removing the driving force for the RFIM behavior in the 
absence of an applied magentic field, that is the param- 
agnetic sublattice. 

Our critical scattering study of CsCo .83Mg . 17 Br3 has 
resolved a longstanding mystery [61 0] as to the ex- 
treme sensitivity of the long range ordered partially- 
paramagnetic Neel state to quenched magnetic vacancies. 
It also contrasts with very recent theoretical work on 
doped quasi-lD Heisenberg antiferromagnets Q. Our 
results show that CsCon.83Mgo.i7Br3 still displays a sym- 
metry breaking near pure CsCoBr'3's Tati^ 28 K, but 
that the system enters a RFIM domain state, and re- 
mains in such a state until sufficiently low tempera- 
tures that the paramagnetic sublattice orders, at approx- 
imately 13 K. 
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